﻿ Stabilitatea sistemelor cu reactie Sistem strict stabil daca – gradul numaratorului lui H(s) (H(z)) § gradul numitorului lui H(s) (H(z)) si/sau – toti polii lui H(s) ∈ semiplanul stang (interiorul cercului unitate pt timp discret) Stabil in sens larg daca polii H(s), sp ∈ axa jω (imaginara) / cercul unitate r=1 pt timp discret polii H(z) zp ∈ Sistemul e stabil numitorul lui H(s) este polinom Hurwitz Q(s) strict/sens larg Hurwitz daca – are coeficienti reali, (strict/nu) pozitivi + – radacini in semiplanul stang (adica au Re{r} 0 In sens larg daca Δk =0 1357 0aaaa 024 6 0aaaa 1350 0aaa 24 60 0naaaΔ= 130 0 0aa 0 a n Ys Hs( Yz Hz ()() )() = (t continuu) = (t discret) 1XzHzGz+ 1XsHsGs+()() () ()() () Ws HsGs= respectiv (Wz HzGz= functia de transfer in bucla deschisa () () ())()() Sistemul in bucla deschisa – strict stabil daca radacinile 1+W(s)=0 ∈ semiplanul stang fara axa jω (partea reala () Hj Gjωω () ,ω∈−∞∞ () inconjoare punctul (-1/k, 0) in sens anti-orar de n ori unde 2nn n= , ic+ in= nr poli semiplanul drept, cn = nr poli de pe axa imaginara (t continuu) 1 jΩ RzHzGz=+ nici un zero sa nu fie in afara cercului unitate r=1 (z=e) ()() () k in= nr zerouri din afara cercului r=1, cn = nr poli de pe cercul r=1 (t discret) 2 Stabilitatea sistemelor cu reactie Probleme Problema 1 Se consideră amplificatorul cu funcţia de transfer Ga () Hs= sa+ tigul în curent continuu al amplificatorului? a) care este câş b) care este constanta de timp a sistemului? te banda amplificatorului ca şi frecvenţa la care modulul răspunsului în c) se defineş frecvenţă al acestui sistem este de 2 ori mai mic decât câş tigul în curent continuu Care este banda amplificatorului considerat? d) Amplificator considerat se plasează într-o buclă de reacţie ca şi în figură: y(t) Gax(t) sa+ K Care este câştigul în curent continuu, constanta de timp, banda amplificatorului cu reacţie? i valoarea lui k pentru care banda sistemului în buclă închisă este e) Determinaţ dublul benzii sistemului în buclă deschisă Rezolvare Ga a) ()Hjω= Câştigul în c c ( 0HG= ) jaω + 1− 1L at− b) ()etσ↔ sa+ t − Dacă () ht e tτσ∼ atunci τ este constanta de timp a sistemului () 1 at− ht Gae tστ=⋅ ⇒= () () a GGa = c) Banda este ()ω== Dar () MBHjHjω 222aω + Ga G aω =⇒= 22T 2aω M+ d) Ys Xs kYs Hs=−⎡⎤ () ()()() ⎣⎦ Ys XsHs kYsHs=− () () () () () 1Ys kHs XsHs+=⎡⎤ () () () () ⎣⎦ 3 Ga Ys Hs ()() sa+ = () bîHs== 1GaXs kHs+ ()() 1k+ sa+ G Câştigul în c c al sistemului cu reacţie este () 0H= bî 1kG+ Ga Ga () bîHs== 1sakGa sa kG++ + + () 1 Constanta de timp τ = 1akG+ () Ga = () bîHω ++ 1jakGω () Ga G ⇒= ++ () + () 121TbîjakGkGω Ga G ⇔= 22 2 21kG+ () ++ () ( ) 1TbîakGω 2 1a 222= 2 ++ + () ( )() 121TbîakG kGω 22 2 22 += ++ ()()() 21 1TbîakG akGω ⇒=+ () 1TbîakGω 1 e) ()akG akωω=⇔+=⇔= 212Tbî T G Problema 2 Determinaţi funcţiile de transfer ale sistemelor din figură a) u(t) x(t) H0(s) H1(s) y(t) G(s) 4 b) H0(z) u[n] v[n] x[n] H1(z) y[n] w[n] G(z) c) x(t) u(t) v(t) w(t) y(t) H2(s) H1(s) G1(s) G2(s) Rezolvare a) ( ) () () 0Us XsH s= ⋅⎡⎤ () () () () () 1Ys Us Ys Gs H s=− ⎣⎦ ⎡⎤ () () () ()() 111Ys GsHs UsHs+= ⎣⎦ ()() 1Ys H s = ()() () 11Us GsH s+ ()() 1Ys H s ⇒= () ()() () 011XsH s GsH s+ ()() () 10Ys H s H s⋅ ⇒= ()() () 11XsGsHs+⋅ b) [ şi () un xn wn=−Vz Uz H z=⋅ ()() ][][] 1 −=− =− () () () ()()()()()()() 1Uz Xz Wz Xz GzVz Xz GzH zUz⇒= ⎡⎤ () () () () 11Uz GzH z Xz+= ⎣⎦ Xz(Vz () ) , dar () Uz=Uz= () () () () 11GzH z+1Hz () () 1HzXz Vz⇒= () () () 11GzH z+ 5 () () () () 0Yz XzH z Vz=+ () () 1HzXz () () 0XzH z=+ () () 11GzH z+ () () 1Yz H z () 0Hz=+ () () () 11XzGzHz+ c) () () () () 2Us Xs G sYs=− () () () 2Vs UsH s= () () () () 1Ws Vs YsG s=− () () () 1Ys WsH s= − ⎡⎤ () () () ()() 11Ys H s Vs YsG s⇒= ⎣⎦ ⎡⎤ () () () () ()() 12 1Ys H s H sUs YsG s=− ⎣⎦ − ⎡⎤ () () () () () () () {} 12 2 1HsH s Xs GsYs YsGs=− ⎣⎦ =⎡⎤ () ()()()()()()()() 122 11 121Ys HsH sGs HsGs HsH sXs++ ⎣⎦ ()() () 12Ys H sH s = ()() () () () () 122 111XsHsHsGsHsGs++ Problema 3 Se dă sistemul descris de ecuaţia diferenţială ′′ ′ ′′′ ′′ ′ 69 68yyyyxxx′′′ +++=++ cu condiţii iniţiale nule Studiaţi stabilitatea acestui system pe baza criteriului Hurwitz Rezolvare 32 3 2 69 6 8sYs sYs sYs Ys sXs sXs sXs+++=+ + () () () ()()()() 32 () 68Yssss++ = 32 691Xss s s++ + () Polinomul de la numitor are toţi coeficienţii pozitivi, deci se poate aplica criteriul Hurwitz = = 01 2 31, 6, 9, 1aaa a== 135610aaa = 3 024190aaaΔ= 130061aa > 160Δ= −= > 254 1 53 0Δ= =⋅⋅+⋅⋅+⋅⋅−⋅⋅−⋅⋅−⋅⋅= =Δ 32691 160 010 090 066 111 53Δ Sistemul considerat este deci stabil 6 Problema 4 Se dă sistemul în buclă închisă din figură x(t) y(t) H(s) k G(s) Determinaţi valorile lui k pentru care acesta este stabil, dacă 1 a () HsGs= () 1s+ 1 HsGs= b () () 1s− 1 HsGs= c () () s⎛⎞ 11s++ () ⎜⎟ 10⎝⎠ 1 d () () 2HsGs= 1s− 1 e () () 2HsGs= 1s+ () 1 f () () 3HsGs= 1s+ () 1 (temă) g () () 4HsGs= 1s+ () 1s+ HsGs= h () () 2 4s− 1 HsGs= i () () 2 22ss++ Rezolvare 11 HsGs= () Punctul a ()ωω⇒= , sp = -1 Hj Gj ()() 1s+ + 1jω Pentru a desena hodograful se folosesc diagramele Bode 1 Hj Gjωω= ()() 2 1ω + 2 ⎛⎞ ⎛⎞ 2ω ⇒ ()20 log 10 log 1 10 log 1Hj Gjωω ω =− + =− + () ⎜⎟ () ⎜⎟ 1⎜⎟ ⎝⎠ ⎝⎠ ω ⎛⎞ arg arg 1 arctg arctgHj Gj jωωω ω =− + =− =− ()()() {} ⎜⎟ 1⎝⎠ 7 Sensul nu e antiorar ω= -1 ω=0 ω=±∞ ω=1 8 Sistemul în buclă deschisă este stabil => hodograful nu trebuie să încercuiască punctul 1⎛⎞ ,0− ⎜⎟ k⎝⎠ 11 00 sau 1 1,0kk⇔− − > ⇒ ∈ − () kk Deci pentru stabilitate este necesar k>-1! 1 Punctul b () HsGs= () 1s− Sistemul în buclă deschisă este instabil cu un pol în semiplanul drept sp=1 1 = Hj Gjωω ()() 2 + 1ω ω⎡⎤ ⎛⎞ arg arg 1 arctg arctgHj Gj jωωωπ πω =− − + =− + =− + ()()() {} ⎜⎟ ⎢⎥ 1− ⎝⎠ ⎣⎦ 111jjωω++ == =− ()() 22Hj Gjωω 11 1jωω ω−−− + 1ω− + 22j=− 11ωω++ Re Im Pentru 0ω > ()Re 0Hj Gjωω 11k−> ⎪ k⎪ Deci sistemul este stabil pentru k>1 1k⇔⇒ > ⎨ 1 − 00k⎪ k⎪ ⎩ 9 Sens antiorar ω=±∞ ω=0 ω=1 10 1 Punctul c () Sistemul în buclă deschisă este stabil HsGs= () s⎛⎞ 11s++ () ⎜⎟ 10⎝⎠ 11 Hj Gjωω = ⇒ ()= Hj Gjωω ()()() 2 jω ⎛⎞ ⎛⎞ 2ω ++ 11jω () ⎜⎟ +⋅+ 11ω ⎜⎟ 10⎝⎠ 10⎝⎠ 22 ⎛⎞⎛ ⎞ ωω ⎛⎞ ⎛ ⎞ =− + − + 20 log 10 log 1 10 log 1Hj Gjωω ()() ⎜⎟⎜ ⎟ ⎜⎟ ⎜ ⎟ ⎜ ⎟ 110⎜⎟ ⎝⎠ ⎝ ⎠ ⎝⎠⎝ ⎠ ω ⎛⎞ =− − arg arctg arctgHj Gjωω ω ()()() {} ⎜⎟ 10⎝⎠ 1⎧ 0− nu trebuie sa fie − Sensul nu este antiorar=>punctul ,0− ⎨ ⎜⎟ 1k⎝⎠ 1⎪ −> k⎪ ⎩ incercuit 11 Sensul nu e antiorar ω=∞ω=0 ω=10 ω=1 1 Punctul d ()1s⇒=± Sistemul în buclă deschisă este instabil HsGs= () 21,2p 1s− 11 ReHj Gj Hj Gjω ωωω==−=; (Im 0Hj Gjωω= ()()()() ) {}{} 22)( 11ωω −− + 0 1 ∞ ωSistemul în buclă deschisă are un pol în semiplanul drept, pentru ca sistemul în buclă închisă să fie Re Hj Gjωω -1 -1/2 0 ()() {} stabil este necesar ca punctul critic să fie încercuit 1 o dată în sens antiorar 101k⇒− k Sensul e antiorar ω=0 ω ω=1 =∞ 12 11 Punctul e ()ωω=⇒ = ()()() 22HsGs H j Gj 11sjω ++ ()() 2 ⎛⎞ 1ω ⎛⎞ Hj Gjωω = () ()() 20 log 20 log 1Hj Gjωω ⇒=−+ ⎜⎟ 2() ⎜⎟ 1ω + 1⎜⎟ ⎝⎠ ⎝⎠ arg 2 arctgHj Gjωωω=− ()()() {} Deoarece sistemul în buclă deschisă este stabil e necesar ca punctul critic să nu fie 11 încercuit 0 sau 1 0 sau 1 0 1kkk⇒− ⇒ > − − kk Sensul nu e antiorar ω=0 ω=10 ω=1 13 11 Punctul f ()ωω=⇒ = ()()() 33HsGs H j Gj 11sjω ++ ()() 23⎡ ⎤ ⎛⎞ 1ω ⎛⎞ ⎢⎥ Hj Gjωω= 20 log 10 log 1Hj Gjωω ⇒=−+ ()() () ⎜⎟ ⎜⎟ 3() ⎜⎟ ⎥ 21⎢ ⎝⎠ ⎝⎠ + 1ω ⎣⎦ ) ( arg 3arctgHj Gjωωω=− ()()() {} 111 1 sau 1 0 sau 8 ,8 1, 0kkk⇒− > − Punctul h () () {} 1,2 122Re 0ppHsGs s s=⇒ 4s− Rezultă că sistemul în buclă deschisă este instabil 11jω+ =⇒ =− Re ,Hj Gj Hj Gjωω ωω ()()()() {} 22 −− + 44ωω ω =− ()() {} ImHj Gjωω2 + 4ω ω ( Re Hj GjωωIm Hj Gjωω ()())() {}{} 0 -1-4 0 1 -1/5 -1/5 ∞ 0 0 Sistemul este stabil pentru 11 0 4 si 0 4kkk− >⇒> 4k Sens antiorar =∞ ω=0 ω ω=1 1 i ()HsGs s j s=⇒=−±⇒ 2kk 0 sau 2kk⇒> >− 2,k⇒∈ −∞ () Sensul nu e antiorar ω=∞ω=0 ω=2 16